Time-dependent Density Functional calculation of e-H scattering 
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Phase shifts for single-channel elastic electron-atom scattering are derived from time-dependent 
density functional theory. The H~ ion is placed in a spherical box, its discrete spectrum found, 
and phase shifts deduced. Exact-exchange yields an excellent approximation to the ground-state 
Kohn-Sham potential, while the adiabatic local density approximation yields good singlet and triplet 
phase shifts. 
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Modern density functional theory (DFT) [J 0, Q 
has proven very succcsful in quantum chemistry and 
solid-state physics. The time-dependent formulation, 
TDDFT [i| , extends this success to excited-state proper- 
ties Thus, excitation energies and oscillator strengths 
of electronic transitions of atoms, molecules, and clusters 
are now routinely calculated via TDDFT within, e.g., the 
adiabatic local density approximation (ALDA) Q. But 
such calculations are almost exclusively for optical re- 
sponse to either weak [H[ or strong Q fields. 

The problem of calculating low-energy elastic electron 
scattering from atoms and molecules is demanding, and 
solving the Schrodinger equation for continuum states of 
polyatomic molecules can be expensive. However, such 
solutions are needed for the emergent field of electron- 
impact chemistry 0] , especially since recent experiments 
show that low energy electrons can cleave DNA [tl, ■ 
Through efficient use of R-matrix theory, calculations 
within static exchange (amounting to scattering from an 
effective one-body potential) have been performed for a 



single DNA base A TDDFT approach could prove 



highly useful here, allowing the incorporation of correla- 
tion effects with little additional cost beyond the original 
scattering calculation. 

With this ultimate goal in mind, we demonstrate a sim- 
ple method for using TDDFT to calculate phase shifts. 
We find the continuum states of the N + 1-electron prob- 
lem, where the target has N electrons. Our method is 
extremely practical in spherical cases, such as atoms. It 
is based on a little-used formula 0, [H, Q [H, EH (ex- 
act for finite-ranged potentials) that relates the phase 
shift of the continuum problem to discrete energies of the 
same potential, but placed inside a box whose edge is be- 
yond the range of the potential. This formula b ypa sses 
many of the complications of our original work [17| , as 
now we need only find bound-bound transition energies, 
where TDDFT has already proven successful. Further- 



more, since our general approach requires that the N + 1- 
electron system be bound, by putting the system in a 
box, our new method can be applied, at least in princi- 
ple, even when the 'ground state' of the TV + 1-electron 
system is only a resonance. 

A vital element in any DFT approach is the accuracy 
of approximate functionals used. In this sense, electron 
scattering from the H atom is a very severe test, since 
H~ (the N + 1-clectron system) is so strongly correlated. 
The underlying ground-state Kohn-Sham (KS) poten- 
tial is crucial to any TDDFT calculation, especially for 
atoms, and is known essentially exactly for H~ [3]. We 
find that exact-exchange, as calculated in an optimized 
effective potential (OEP) code, yields very accurate KS 
phase-shifts, i.e., very close to those of the known ex- 
act KS potential. Next, we show that the ALDA, the 
workhorse of TDDFT, yields very good shifts for both 
singlet and triplet (TD-spin-DFT) scattering. Thus, we 
demonstrate that a simple formalism allows scattering 
calculations from TDDFT; that modern approximations 
yield sufficiently accurate ground-state potentials; and 
that standard TDDFT approximations are sufficiently 
accurate. We perform the first such calculation on the 
prototype target, the H atom. 

We begin with some exact observations about scat- 
tering from a potential. Consider a spherical potential 
that has a finite range, i.e., v(r) = beyond some radius 
R c . Now imagine inserting a hard wall at any Rb > R c , 
not necessarily far beyond R c , and solving for the bound 
states. Any such solution is in fact a solution to the origi- 
nal scattering problem that happens to have a node right 
at Rb- Study of the wavefunction between R c and Rb to 
identify the phase shift yields: 



tan(<5 iQ ) = -ji(k a Rb)/m(k a Rb) 



(1) 



where ji and r]i are the two free-space solutions to the 
radial Schrodinger equation, i.e., the spherical Bcsscl and 



2 



e-H phase shift 



■ Accurate data 
23 a.ii. 
i 100 a.u. 




H' potentials 




> -0.1 - 



FIG. 1: Accurate quantum chemical singlet and triplet s- 
phase shifts , together with the KS values, calculated with 
a wall at 23 a.u. and at 100 a.u.. 



von Neumann functions, k a = \/2E a , and E a is the cv-th 
eigenenergy. For s-wave scattering, Eq. (p} reduces to: 



5a — k a Rb 



(1 = 0) 



(2) 



The phase shifts are only determined modulo 7r, but we 
have added an, the free-particle value of k a Rb, so that 
all shifts are relative to 0. For any given I and Rb, this 
method yields the phase shift at a discrete set of energies. 
For a fixed potential, starting from any Rb value, i?n, one 
can continuously increase Rb to about 2Rq and generate <5 
at all energies above a minimum £" m i n = ^[(7r — 8 n )/Rb] 2 , 
where Rb is the largest box used. Usually i?-matrix the- 
ory is more convenient, as it does not require the wave- 
function to have a node at the box radius, and so all 
energies can be found with just one value of Rb- But it 
relies on knowing the logarithmic derivative of the wave- 
function at Rb, which is not available in TDDFT. 

To illustrate the method, and show how useful the ex- 
act ground state KS potential is, in Fig. [T]we plot ac- 
curate quantum calculations for both singlet and triplet 
elastic scattering from hydrogen (l9| . We also plot the 
result of potential scattering from the exact ground state 
KS potential of H~. This was found by Umrigar et 
al. [18|, from an extremely accurate quantum Monte- 
Carlo calculation for the ground state of H~ , calculating 
the density, and finding v s (r) by inverting the KS equa- 
tion. We obtained the positive orbital energies (necessary 
to evaluate Eq. ©) from a well-established fully numer- 
ical spherical DFT code, which includes the optimized 
effective potential method (OEP) and has been supple- 
mented by the option to insert a hard-wall at a distance 
Rb from the origin [20j|. 

The KS phase shift fits between the two curves, just as 
the pure KS orbital energy differences lie between singlet 
and triplet excitations for He (21, 22 1. The calculations 
at two values of Rb demonstrate the results are indepen- 



FIG. 2: The exact, exact-exchange, and LDA KS-potentials 
for ET. 
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FIG. 3: The s-wave phase shifts for the exact and exact- 
exchange KS potentials for H~. 



dent of Rb- We choose the wall far from the origin to 
ensure the self-consistent ground-state results are not af- 
fected by its position and to approach zero energy, but we 
emphasize the fact that Eq. is exact for any Rb > R c - 
Obviously, the decay of the ground-state KS potential 
is crucial to the accuracy of this method, and in any prac- 
tical application, the exact KS potential is unavailable. 
Therefore we study the behavior of two approximate po- 
tentials, exact- exchange (OEP) and the local density ap- 
proximation (LDA) . In Fig. [2J we plot both the exact and 
approximate KS potentials for H~ . The LDA potential is 
far too shallow, a well-known failing of most commonly- 
used approximations to ground-state DFT. The true self- 
consistent LDA potential does not support any bound 
states, so to obtain the potential we put the system in 
a large box, forcing the states to be bound [23(. Thus, 
the LDA potential is utterly unsuitable for this type of 
calculation. On the other hand, the exact-exchange po- 
tential decays correctly as r — > oo, missing only the small 
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positive correlation potential for small r. Many modern 
i?-matrix based methods start with the nuclear potential 
and the pure electrostatic (i.e., Hartree) potential, and 
then add the LDA exchange potential from DFT, i.e., the 
Slater contribution that decays exponentially, as n 1 / 3 (r). 
Since this potential misses the correct asymptotic behav- 
ior, fxcJz) —* —l/ r , a 'polarization' potential must be 
added [24|. Our KS potentials, either exact or exact- 
exchange, already have the correct asymptotic behavior, 
i.e., they contain the polarization potential. Without this 
feature, our KS potentials would have the wrong asymp- 
totic behavior, and would not be long-ranged for neutral 
atoms. In Fig. [3j we plot the scattering from the ex- 
act and exact-exchange potentials, demonstrating that 
exact-exchange, as is now available in many codes [iif . 
is perfectly adequate for this purpose. 

To go even further, e.g. to account for singlet-triplet 
splitting, we must use TDDFT. Within the formalism of 
TDDFT within linear response we can, in principle, ob- 
tain the true singlet and triplet excitation energies, and 
thus the phase shifts. We label all single-particle excita- 
tions from the ground to unoccupied excited states via 
q = (i,a), where i implies occupied, a implies unoccu- 
pied, and define $ 9CT (r) = 0* CT (r)</) a(T (r), where a is a 
spin index and 4>j(r) is an eigenstate of the ground state 
i> s (r). Casida |26[ cast the TDDFT response equations 
as a eigenvalue equation 
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FIG. 4: Singlet and triplet TDDFT curves from an SPA and 
full ALDA calculation, together with the KS values and ac- 
curate quantum chemical data from Ref. The ground 
state KS potential is exact-exchange. The wall location in all 
calculations is at 100 a.u. 



TABLE I: TDDFT s-wave scattering lengths. 
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Triplet a 


Accurate data" 


5.97 


1.77 


ALDA SPA 
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1.8 


ALDA 


5.6 
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"Accurate variational calculations from fl9ll 

''This is the value as obtained from our tangent approximation as 
explained in the text 



where 

Qqo-q'cr' (w) = LO^ a 5 qq i 8 aa i 

+ 2^LU qa oj q , a ,(qa\f^' c (uj)\q'a'). (4) 

and (qa\f^' c (us)\q'a') is the matrix element of the 
Hartree- XC kernel in the set of functions $ g(T (r) . We also 
defined us qa = €i a — e aa , where ej CT is the KS orbital energy 
of state i with spin a. The XC kernel is the functional 
derivative of the XC potential in TDDFT [!, Q and we 
assume the frequency independent ALDA kernel in the 
following. The f2-matrix can be split in separate singlet 
and triplet tl- matrices (26j . Solving Eq. §3§ therefore 
yields predictions of both singlet and triplet transition 
frequencies, us. In order to perform these calculations we 
have added subroutines to evaluate the matrix elements 
needed for a TDDFT calculation. Since the system stud- 
ied is small, we exactly diagonalize the f2-matrix. 

In Fig. QJ we show the results obtained from a 
TDDFT ALDA calculation, but using the OEP ground- 
state potential. Apart from the full results we also 
show results obtained with the single-pole approximation 
(SPA), which ignores the off-diagonal matrix elements in 
Eq. (0} [23] . The SPA is analogous to the distorted- wave 
Born approximation used in earlier work, which worked 



well for electron scattering from He + [l7| , but fails badly 
for H. In Fig. [4] the SPA is indeed a poor approxima- 
tion to the full singlet curve, especially at low energies. 
For the triplets, on the other hand, we obtain excellent 
results that are on top of the reference values. 

If we now look at the full calculation we see that in- 
cluding the off-diagonal matrix elements considerably im- 
proves the singlet values giving results very close to the 
reference data. For the triplet the results only change for 
smaller energies where the values are too big and there 
is a small "bump" close to E = 0. We believe this bump 
to be unphysical, due to coupling among transitions be- 
ing treated incorrectly by our approximate XC-kernel, as 
E — > 0. This suspicion is reinforced by the fact that in 
this region, the full ALDA triplet results depend on the 
position of the wall, and so cannot be trusted. However, 
similar effects were found with other common kernels, 
such as exact exchange, so we believe some delicate be- 
havior of the XC kernel is required to avoid this artifact. 

To quantify results for low energies, the scattering 
length is defined by the effective range expansion, 



k 2l+1 cot 5i(k) = - — + -r el k 2 
y 'k^o at 2 



0(fc 4 ), (5) 

where ai is the scattering length and r e ; the effective 
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range. Since we have no wave function, to extract a; we 
must fit our data to the above expression to obtain the 
scattering lengths. We give a rough estimate of the ex- 
pected TDDFT scattering lengths in Table [fl by fitting 
our results for small k. As a reference, the KS scatter- 
ing lengths are 4.7 for the exact potential and 4.2 for 
the exact-exchange potential. We report the scattering 
lengths we obtained from our phase shifts, or in the case 
of the full triplet ALDA calculation, from fitting a tan- 
gent to the curve, required to pass through ir. In the 
triplet case, the value obtained from the fit agrees well 
with the reference value, as does the SPA result. Thus 
either method yields accurate results as k — > 0. 

Scattering from neutrals is very different from scatter- 
ing from positive ions. In the former, the N + 1-electron 
system has a s/iori-ranged potential, and so a finite cross- 
section, but in the latter, the KS potential is long-ranged, 
i.e., it decays as —1/r for large r, and the cross-section 
diverges. The phase-shift is then defined relative to pure 
Coulomb scattering. Our general approach still applies, 
but Eq. (fT]) must be modified. If a potential deviates 
from —1/r only for r < R c , 



sions with Mike Morrison, Neepa Maitra, Chris Greene, 
Peter Lambropolus and Vladimir Mandclshtam. 



tan(<5 Q ) = -Fi{k a Rc)/Gi(k a Rc) 



(6) 



where Fi and Gi are the Coulomb scattering solu- 
tions [28|. We will report results for positive ions in 
future publications. 

While the ALDA functional uses only input from the 
uniform electron gas, our results show that it gives ac- 
curate results for electron-scattering from a system that 
could not be further away from a homogeneous gas, the 
hydrogen atom. These results encourage us to continue 
work along these lines. We will calculate other Z-values, 
different approximate ground-state potentials, different 
XC-kernels, other atoms, and ions, to gain experience 
in the reliability of TDDFT calculations. But we finish 
by considering some obstacles in applying our method to 
scattering from large molecules. We first note that, by 
converting the problem to one of discrete transitions, one 
needs only modify an existing electronic structure code 
by placing a hard wall around it, rather than use a scat- 
tering code. However, our formula is only exact if the 
wavefunction has a node on the hard-wall surface, which 
would only be true state-by-state for a non-spherical sys- 
tem. Much better is to use a large sphere, so that the 
formula is approximately true. We must also address 
the multichannel case. We intend applying our method 
to electron scattering from Be + next, which is a well- 
studied scattering example (29[, and for which the exact 
ground-state KS potential is known [l8j |. 
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